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Abstract 

We discuss the gravitational lensing of neutrinos by astrophysical 
objects. Unlike photons, neutrinos can cross a stellar core; as a result, 
the lens quality improves. We also estimate the depletion of the neu- 
trino flux after crossing a massive object and the signal amplification 
expected. While Uranians alone would benefit from this effect in the 
Sun, similar effects could be considered for binary systems. 
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1 Introduction 



In this letter we investigate the possibility of neutrino gravitational lensing 
by astrophysical objects (stars, galaxies, or rather galactic halos). Unlike 
photons, neutrinos can cross even a stellar core. This results in a much better 
focalization due to the improvement in the lens quality. However, because of 
the extreme difficulty of detecting neutrinos and the poor angular resolution 
of "neutrino telescopes", one can only expect signal intensification rather 
than a spectacular lensing pattern like those of photons. The possibility 
of neutrino focusing by the Sun's core was previously debated in Ref. 0. 
Neutrino lensing has also been used as a possible explanation for the time 
difference between neutrinos from supernova SN1987a ||. Here, we pursue 
an exhaustive and general analysis of neutrino lensing involving calculation 
of the neutrino deflection with different density profiles, a comprehensive 
discussion of neutrino absorption and signal magnification, and a description 
of possible applications. The letter is organized as follows. The gravitational 
deflection of neutrinos is calculated in Section ||. In Sect. |3], we study the 
neutrino interactions in a stellar medium and the depletion of the neutrino 
flux after crossing a massive object. Signal enhancement is estimated in 
Sect. Finally, in Sect. [5], we present some practical examples of our analysis. 
It is clearly shown that the Earth-Sun distance is too small for a sizable 
lensing effect to take place, although an observatory on Uranus would notice 
the enhancement of distant neutrino sources whenever they are aligned with 
the Sun. The case of galaxies, or rather their halos, is also contemplated. We 
finally consider binary systems, which are the most promising case. A more 
extensive and detailed description of our investigation, including complete 
calculations of the gravitational lensing effect, neutrino absorption and signal 
enhancement (taking into account the geometry, the amplification factor, and 
certain approximations of the lensing phenomena) can be found in Ref. Pfl. 



2 Gravitational deflection of neutrinos 

In this section we study the deflection of neutrinos from straight-line motion 
as they pass through a gravitational field produced by a compact object 
of mass M and physical radius R. We will distinguish two cases: when 
the neutrino flux passes far away from the object (OUTside solution), a 
situation equivalent to the gravitational lensing of photons, and when the 
neutrino flux passes through the object (INside solution). In the latter case, 
we consider three specific cases depending on the compact object density 
profile: constant density (a profile fit for planetary objects), Gaussian density 
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distribution (suitable for starsQ) and Lorentzian density distribution (which 
could be associated with a galactic halo|). 

For the OUTside solution, calculating the trajectory of a massless neu- 
trino (or with a mass very small compared to its energy) in the Schwarzschild 
metric under the assumption that M/r is everywhere small along the trajec- 
tory ||, U gives as a result for the net deflection angle A0out = 4M/&. For 
the case of a neutrino flux passing through the object, one must first, in 
order to study the neutrino trajectory, look for the form of the space-time in 
the region inside the object. Here we restrict ourselves to the case of static 
spherically symmetric space-times and static perfect fluids. In the region 
inside the object, exact solutions to the relativistic equations are very hard 
to solve analytically for a given equation of state M. One interesting exact 
solution is the Schwarzschild const ant- density inside solution, which we use 
here as a typical density profile for planetary objects. The net deflection as 
a function of the impact parameter b is then 
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where the outside solution is also included for completeness. In Fig. [I], we plot 
the normalized net deflection A0/A0(6 = R) as a function of the normalized 
impact parameter b/R for the constant distribution density and compare it 
with other specific density profiles considered in the analysis. Such a normal- 
ization allows for a clear comparison of different profiles and is independent 
of the mass and physical radius of the compact object. 

We refer the reader to Ref. || for a detailed description of the calculations 
done in this section. 

Next, we analyze the solutions for the Gaussian and Lorentzian densities 
distribution. The Gaussian profile is a convenient approximation of the mass 
distribution in stars while the Lorentzian profile is valid for galactic halos. 
In both cases, it is possible to neglect the pressure with respect to the mass 
density, p <^ p (see Ref. [|J for the so-called Newtonian stars), and thus 
4wr 3 p <C m. For the case of a Gaussian density profile p(r) = po e_r ^ r °> the 

2 The density profile of stars is not exactly Gaussian but we use it here in order to obtain 
simple analytical results. Such a description should be considered as a good approximation 
to the real case. 

3 The Lorentzian profile behaves as 1/r 2 for large r, in agreement with velocity disper- 
sion curves for galaxies and clusters. 
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Figure 1: Normalized net deflection A0/A0(5 = R) as a function of the 
normalized impact parameter b/R for a constant density profile (solid line), 
a Gaussian density profile (dashed line) and a Lorentzian density profile 
(dotted line). In the last two cases r is taken to be r = 0.2R. 
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where the error function is defined as erf (z) = ^= Jq dte * 2 . Fig. [T] shows 
the exact result in Eq. (0) with the parameter r fixed to r = 0.2i?. As 
it is seen from Fig. the maximal net deflection occurs at b ~ r and is 
A0(6 ~ ro) — 3.2A0(6 = i?), showing that the lensing effect inside the star 
is bigger than the outside effect (except for b < 0.04_R). 



For the case of a Lorentzian density profile p{r) 
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Figure 2: Normalized focal length f/f(b 
impact parameter 6/i?. 
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flection is 
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In Fig. we also plot the exact result in Eq. (|3|) again with r = 0.2R. 

In summary, we conclude that only for the case of stars (Gaussian profile) 
is the lensing effect inside the star substantially amplified with respect to the 
outside effect (the inside effect should be compared with the effect at b = R) . 
For galactic halos, the maximal inside net deflection is slightly bigger than 
at b = R, while for an object of constant density the inside lensing effect is 
always smaller than at b = R. 

Finally, we plot the focal length / as a function of the impact parameter 
b. The effective focal length is defined as the distance at which the lens 
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focuses the signal 

'<»> = a* ■ (4) 

A perfect lens would correspond to a constant focal length (independence 
of b). The neutrino flux would then be focalized in a single point with 
consequent signal intensification. In terms of the net deflection, a "good 
lens" requires that A0 increases with the impact parameter b. In Fig. |2|, 
the (normalized) focal lengths for the three different profiles are shown. We 
postpone to Sections. |] and |5] our comments about the quality of the different 
gravitational lenses considered here. 



3 Neutrino absorption 

As long as neutrinos are able to travel across a massive object, one has 
to consider their interactions with the matter inside. These interactions, 
dominated by scattering on nucleons || , will indeed reduce the neutrino flux 
and thus the efficiency of signal amplification due to lensing. Charged-current 
reactions convert the neutrinos into charged leptons that will later decay 
or be absorbed by matter while neutral-current reactions deviate them by 
large angles compared to the deflection angle. In both cases, the interacting 
neutrinos will not contribute to signal enhancement. 

The probability of transmission for a neutrino crossing at a distance b 
from the center of the object isQ 

P T (b, E u ) = exp (-2 a uN (E u ) j* N N (r) jj^ j , (5) 

where a uN (E u ) is the neutrino- nucleon cross section as a function of the 
laboratory neutrino energy and Nn is the number density of scatterers for 
the case of nucleons. In order to compute the probability of transmission, one 
has to collect information not only on the scattering cross section but also 
on the mass density profile and the composition of the object (the latter is 
needed in order to relate number and mass densities). A detailed calculation 
of Eq. (D is found in Ref. §. 

Let us now briefly discuss the results and consequences of Eq. (H) for 
the two relevant physical cases studied in Sect. i.e. neutrinos crossing a 
star or a galactic halo. For the case of a star, Eq. (|) is calculated taking 

4 We assume for simplicity an object with spherical symmetry. 
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Figure 3: Neutrino probability of transmission PT(b,E u ) as a function of 
b/R for neutrinos crossing the Sun with energies E v = 10 GeV, 10 2 GeV, 
10 3 GeV, 10 4 GeV, 10 5 GeV and 10 6 GeV (solid lines from left to right). 
The Gaussian density profile with r$ = 0.2R & is included for completeness 
(dashed line). 

the Sun as example^]. The probability of transmission of neutrinos^ across 
the Sun is then shown in Fig. ^| as a function of b/R for different incoming 
neutrino energies. From Fig. |3|, it comes out that E v > 1 TeV neutrinos are 
completely absorbed or scattered away in the core of the star. Even at 100 
GeV, matter interactions are frequent enough to significantly attenuate the 
flux leaving the star. Consequently, the region inside the star where lensing 
is more efficient (it was shown in Sect. § that the star's core acts as a "good 
lens") is ruled out as a lens for neutrinos with E v > 100 GeV due to flux 
attenuation. For E v < 100 GeV neutrinos, however, matter interactions have 
practically no incidence on the outgoing flux which is then recovered after 
lensing by the star. 

Neutrinos passing through a galaxy may interact either with its visible 
matter or with the surrounding halo of massive relic neutrinos. For the 
former, a rough estimate of the average density of stars in a galaxy gives 
1 pc~ 3 , which results in a negligible probability for a neutrino to encounter 

5 Neutrino flux attenuation in the Sun has also been considered in Ref. jjj. 

6 A similar figure is obtained for antineutrinos. In this case, however, there also exists a 
contribution from v e e scattering that dominates over the vN contribution at the Glashow 
resonance, E^ s — M^/2m e w 6.3 x 10 6 GeV, but it is negligible at all other energies 
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a star during its passage through the galaxy, even if it traverses the whole 
disk and the bulge. For the latter, only interactions with ultrahigh energy 
neutrinos are significant |3J and thus are not taken into account in our present 
framework. We conclude that the passage of neutrinos through a galaxy will 
not decrease their flux and hence not influence the lensing effect at all. 



4 Signal enhancement 

The difficulty of detecting neutrinos together with the poor angular resolution 
of "neutrino telescopes" leaves signal enhancement as the most likely (if not 
only) signature for neutrino lensing. Here we discuss in some detail this 
signal enhancement for different neutrino sources and lenses. 

The total amplification (or magnification) of the signal intensity fi for a 
point source is given by B iKl 



pdp 



(6) 



where Ai and Aq are the surfaces of a certain image (in general there is more 
than one) and the source (both referred to the lens plane), 6i are the angles 
of the source images with respect to the observer- lens line of sight, and (3 is 
the angle to the actual source position. For a specific image with angle 9, 

P = 9-^A<P, (7) 

-^so 

A(f) being the deflection angle and D\g(D so ) the distance between the source 
and the lens (observer) (see Fig. f|). 

To observe neutrino lensing would require a huge amplification. As de- 
duced from Eq. this occurs only when the source, the lens, and the 
observer are in perfect alignment, f3 = 0, or nearly aligned, (3 -C 9 S , where 
9 S is the angular size of the source. When alignment occurs, the finite size 
of the source cannot be neglected and the amplification is then calculated as 
i 

^ = -h t m^m- (8) 

ir9g Jo 

In Eq. (H) the detector is taken as pointlike. To consider its finite size[] 
would be relevant only for R d 3> R s %^, where R^Rs) is the size of the 
detector (source), that is for the case of very small sources or very nearby 

7 The case of a finite sized detector is described in Ref. H . 



7 



p 



* ► 

D LS D lo 

Figure 4: The geometry of a lensing event. The lens is located at a distance 
D\ Q from the observer and the source at a distance D so (the distance between 
the lens and the source is D^). The angular separation between the lens and 
the source is (3 and the position of the image is at an angle 9. The deflection 
angle is a (notice that we use A<f> for the deflection angle along the main 
text). 



lenses. These finite size effects should be taken into account for the Sun in 
case it could act as a lens. 

Due to the complicated dependence of the net deflection A0 on the impact 
parameter b (or, in terms of 9, trough the relation 9 = b/D\ Q ) for the different 
density profiles discussed in Sect. |], the amplification //(/?) as well as the 
function 9(/3) are only known implicitly. The only exception is of course the 
OUTcase, where the deflection angle is simply A<p = 4M/b. For this reason, 
we calculate here the signal magnification for two instructive approximations 
of A0 in the INcase. The OUTcase is also shown for the sake of completeness. 
The results are the following: 

A0 = 4M5, n(p) = *fr, * = (9) 
A0 = 4M 7 6, ^(/5) = (^) 2 , ;, S =(W) 2 , 

where 5 and 7 are constants with respect to b, 9q is the value of 9 for perfect 
alignment (9 is nothing but the Einstein's radius, Re = 2-y // MJ ^'° D -^, in an- 
gular units), and #disk is the angle for which the linear approximation in A<p 
holds. All the previous results are valid provided /3, 9 S <C 9 , #disk- The sec- 
ond column presents the magnification of a point source while the third one 
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is for a finite sized source and for perfect alignment. The sum over multiple 
images is included. 

As stated in Eq. (^), the magnification for a constant deflection is just 
twice that of the OUTcase. For a linearly increasing deflection the resulting 
magnification is instead the square of the OUTcase. The former case applies 
rather well to narrow Lorentzian profiles (small width, r$ <C R, typical of 
galactic halos) while the latter fits correctly the central regions of the con- 
stant and Gaussian density profiles. This last case corresponds to the best 
amplification we can hope for since the signal crossing the central region is 
focused. This focalization is expected as soon as the focal length is almost 
constant in the central region, as seen in Fig. ^]. Accordingly, the central 
region of planetary^ and stellar objects acts as a "good lens" in the optical 
sense. Unfortunately, planetary objects are generally too small and light to 
be of practical use. Stars, however, can provide a huge signal amplification 
as well, as they also have a significant "good lens" region. We will now 
investigate in some detail the possible applications. 

5 Applications 

As we are looking for huge amplifications, we consider the case of perfect 
alignment. f(b) is the effective focal length for impact parameter b. Focal- 
ization on the observer will take place if f(b) = f = (-^- + -jj-") • A good 
lens corresponds to a constant f(b). 

The Sun: We use as announced a Gaussian density profile as a reasonable 
approximation of matter distribution in the Sun. Unfortunately, when 
f(b) is compared with / ~ D\ = -D©_® = 1 au (as D\ s ^> D\ is 
assumed), it appears that there is no intersection between the two 
curves. This means that the Sun cannot focus on the Earth a neutrino 
beam coming from a far source. The required intersection happens 
nevertheless at distances around 20 au, meaning that Uranians could 
perform neutrino lensing experiments using the Sun as a lens (see also 
Ref. HI). For them, any neutrino source would be amplified in turn as 
the Sun sweeps in front of it! It is easy to check that Jupiter cannot 
replace the Sun as a useful lens for us, as its mass is about 1O _3 M . 

If one is not located at the focal point there is yet some amplification 

8 The constant profile focuses a beam the most efficiently, as the focal length remains 
approximately constant over a wide range of b/R. 
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provided by the expression^ 



u Gauss = ( ffl i no) 

Eq. (|T0D shows that in order to achieve an amplification factor greater 
than 2, D\ must be in the range (0.3-1.7)/(6). Outside this interval 
we can consider the amplification negligible. 

We conclude that the Sun is not a neutrino signal amplifier for experi- 
ments on Earth. 

Stars, the lighthouse possibility: Stars are notorious lens candidates. 
Each alignment of a neutrino source with a star will provide a lensing 
event, i.e. there exists an intersection between f(b) and /. However, as 
Ac is in this case at least of 1 pc (which is indeed the distance to the 
closest stars), f(b) will always cross / well outside the starQ. The case 
is thus similar to that of photons, namely the lens focuses a thin ring 
whose radius is the Einstein radius. According to Eq. (^), a typical 
magnification factor for stars would be (assuming D\ s 3> D\ ) 



/i s OUT ^4x!0 4 



R S J v 10kpc/ \ \M & ) \ D lo 



M \ /lOOpc 



where _R S is the physical radius of the source. Even if huge amplifica- 
tions are possible, we have verified that for expected neutrino sources 
this signal enhancement is insufficient to allow the detection, the reason 
being that the magnification cannot compensate the 1/D^ Q geometrical 
suppression of the signal. 



We shall now discuss the case of binary systems ||10|| , where one star 
acts as a source and the other, very close to the first one, acts as a lens. 
Binary systems satisfy D\ s <C D\ Q and / ~ D\ s . The most interesting 
behaviour occurs when the source is located at the effective focal length 
dictated by the central region of the lens; the image is then a full disk. 
This coincidence happens for 



1 , (12) 




9 We retain the notation f(b) in Eq. ( |l0| ) so as not to be confused with / although in 
this case f(b) is a constant. 
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A simple calculation taking a Sun- like star as example and D\ Q — 0(1 pc) gives b 



^JAMqD\ ~ 20i? Q for the impact parameter. 
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where f(b) = Tq/AM in the central region of a Gaussian profile and the 
width is set at r = 0.2R. This condition can clearly be fulfilled. An 
upper estimate^ of the total magnification is given by (see Ref. M for 
details) 

where ^disk is the angular size of the disk image. The source should 
then be at least one order of magnitude smaller in radius than the lens 
in order to provide a large magnification. Binary stars will seldom meet 
all the above requisites but for more exotic systems, with a compact 
and intense neutrino source, this situation is really promising. In some 
way, the large companion acts as the lens focusing a lighthouse beam. 

Galaxies: Lastly, we shall discuss another class of lens candidates, namely 
galaxies or rather galactic halos. In Fig. [|, it is shown that the maxi- 
mum deflection angle for a Lorentzian density profile is close to AM/ R 
and, more precisely, is always smaller than ir/2 -4M/R (independently 
of the value taken for tq). So, in order to focus on Earth, an effective 
focal length of 

^-0O P eg)(^) 2 (^), (14 , 

is needed. Therefore, since the radius of the Universe is about 5 Gpc, 
the neutrinos have to pass near the center of the galaxy, often through 
its visible part. Unlike photons, which are absorbed by dust clouds, 
neutrinos survive in their travel through the galaxy because the proba- 
bility of meeting a star is tiny. For calculating the signal enhancement 
we use the Lorentzian density profile. In this case, however, the precise 
value of the width r is not known, although it has to be small enough 
to explain the velocity dispersion curves inside galaxies. Assuming r 
is small we calculate the magnification for b > r$ (this time the image 
is a ring): 

As in the case of binary systems, a small neutrino source is required 
for non- negligible signal enhancement. 

11 The focal length is not exactly constant over the whole central region of the star but 
suffers from hyperbolic aberrations that indeed reduce the total magnification M. 
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Experimental prospects: Neutrino telescopes are putting limits on neu- 
trino fluxes from point sources. Amanda quotes an upper limit on 
Earth of $„ S upper = 10~ 7 cm -2 s _1 , assuming an E~ 2 spectrum, 



for E v > 10 GeV and a declination larger than 30 degrees |Tl], [12 
Active galactic nuclei (AGNs), gamma-ray bursts (GRBs), supernova 
remnants (SNRs), emissions from accretion disks in binary systems, 
etc. are candidates for such sources (for a review, see Refs. |13[ H). 



Gravitational lensing increases the sensitivity to these sources but acts 
only fortuitously. 

A permanent amplification /i (we mean here that, relative to the source, 
the angular speed of the lens is negligible at our time scales) increases 
the expected number of events and therefore the sensitivity by the same 
factor fi. However, to confirm that gravitational lensing is acting, the 
multiplicity of the image must be observed; it thus requires a photonic 
counterpart. 

Non-permanent lensing events will only be observed if enough interac- 
tions are tracked in the detector, i.e. if J* lens ev ji(t) t dt & u > T$ upper , 
where T is the total observation time and $ up per refers to the upper 
limit on point source flux reached by the experiment. The amplification 
pattern is in principle sufficient to demonstrate the gravitational lens- 
ing but requires enough statistics to establish the shape of the event. 
A coincidence with the photonic counterpart can again confirm the 
lensing effect. 

Finally, binary systems allow for periodic amplifications. The average 
sensitivity is increased by a factor f Tev ji(t) t dt and will in most cases 
improve only slightly (as the alignment time is always much shorter 
than a revolution). If the detector is not sensitive enough to match the 
unamplified source, it may collect the amplified signal periodically. It is 
worth pointing out that the neutrino signal is amplified while the source 
is hidden; the photon and neutrino fluxes have thus anticorrelated time 
dependences. It should be noted however that small sources are needed 
for efficient amplification. 



6 Conclusions 

In this letter we have examined neutrino gravitational lensing by astrophys- 
ical objects. Unlike photons, neutrinos can cross a stellar core. We have 
calculated the deflection angle and the effective focal length for different 
density profiles. The neutrino absorption is also discussed and signal ampli- 
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fication is estimated for the relevant cases. Our formalism is then applied to 
stars and galaxies (or galactic halos). 

Our analysis shows that neutrinos passing through the central region of a 
star are deflected at an angle increasing linearly with the impact parameter. 
This results in a "good lens" that can focus neutrinos into a small region. 
The phenomenon operates up to large energies (about 300 GeV for the Sun) 
before the star becomes opaque to neutrinos. Unfortunately, focalization by 
the Sun occurs only at the distance of Uranus, so that amplification on Earth 
is modest. A more promising case is that of binary systems, where a large 
object would focus and "beam" neutrinos originating from a smaller source. 

Acknowledgements 

Work partly supported by IISN Belgium and by the Communaute Frangaise 
de Belgique (Direction de la Recherche Scientifique programme ARC). V. Van 
Elewyck and D. Monderen benefit from a FRIA grant. 

References 



[i] 

[2] 
[3 



[4 

[5 
[6 

[7 



J. L. Gerver, Phys. Lett. A 127 (1988) 301. 

J. D. Barrow and K. Subramanian, Nature 327 (1987) 375. 

R. Escribano, J. M. Frere, D. Monderen and V. Van Elewyck, |hcp- 



ph/9910510 



B. F. Schutz, A first course in general relativity, Cambridge University 
Press 1985. 

S. Weinberg, Gravitation and cosmology, John Wiley & Sons 1972. 

R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Phys. Rev. D 58 
(1998) 093009 ||hep-ph/9807264 |. 



G. Ingelman and M. Thunman, Phys. Rev. D 54 (1996) 4385 fhep- 
ph/9604288| 



E. Roulet, Phys. Rev. D 47 (1993) 5247. 

[9] P. Schneider, J. Ehlers and E. E. Falco, Gravitational lenses, A&A Li- 
brary, Springer 1999. 



13 



[10] E. Roulet and S. Mollerach, Phys. Rept. 279 (1997) 67 ptro= 
| ph/9603TT9l . 



[11] E. Andres et al. [AMANDA Collaboration], Nucl. Phys. Proc. Suppl. 91 
(2001) 423 ||astro-ph/000924l] . 

[12] J. H. Kim [AMANDA Collaboration], 26th International Cosmic 
Ray Conference (ICRC 99), Salt Lake City, UT, 17-25 Aug 1999, 
[http : / /krusty . physics . utah . edu/~icrcl999/proceedings . html] . 

[13] T. K. Gaisser, F. Halzen and T. Stanev, Phys. Rept. 258 (1995) 173 
[Erratum-ibid. 271 (1995) 355] ||hep-ph/ 94 10381 - 

[14] R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Astropart. Phys. 5 
(1996) 81 ||hep-ph/9512364| . 



14 



